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Abstract 

We report the hrst experimental demonstration of combined spatial and temporal control of light 
trajectories through opaque media. This control is achieved by solely manipulating spatial degrees 
of freedom of the incident wavefront. As an application, we demonstrate that the present approach 
is capable to form bandwidth-limited ultrashort pulses from the otherwise randomly transmitted 
light with a controllable interaction time of the pulses with the medium. Our approach provides 
a new tool for fundamental studies of light propagation in complex media and has potential for 
applications for coherent control, sensing and imaging in nano- and biophotonics. 



Concentrating light in time and space is critical for many applications of laser light. 
Broad-band mode-locked lasers provide the required ultrashort light pulses for multipho- 
ton imaging , nanosurgery [3j, microstructuring [4J, ultrafast spectroscopy and 
coherent control of molecular dynamics or of nanooptical fields [z-9]. Multiple random scat- 
tering in complex media severely limits the performance of these methods, but often is an 
unavoidable nuisance in many systems of interest, such as biological tissue or nanophotonic 
St _ Q. Spatial r a„do m sca tt e ring — d ist0rt s a propagating wave front, 
creating the well-known speckle interference pattern In the time domain, ultrashort 

pulses are strongly distorted and widely stretched due to the broad path length distribution 
in multiple scattering media [ijj]. These temporal and spatial distortions are not separable 

There is a strong interest in improving applications of ultrashort laser pulses in com- 
plex scattering media. Phase conjugation has been applied to spatially focus light from a 
short-pulse laser source through a thin scattering layer [14]. Similarly, phase conjugation 
is applied to correct distortions of the ballistic wave front to improve the resolution of two 
photon microscopy [l5| . Coherent control of two-photon excitation through scattering bio- 



logical tissue has been demonstrated 



16]. Those experiments share the common limitation 



that the control is limited only to those photons that take the shortest paths through the 
disordered media and arrive a the target volume without being multiply scattered. 

Recently it was demonstrated that random scattering can actually be beneficial rather 
than detrimental for the performance of optical systems. Applying a shaped wave front of 
monochromatic light to a str ong ly scattering medium, Vellekoop et al. achieved spatial con- 
trol over the scattered light 17] . In fact, the insertion of an opaque sample after a lens has 



allowed focusing beyond the diffraction limit of the lens 



18] . These findings have opened new 



possibilities for imaging in optically thick biological matter 19] and allow trapping particles 
through turbid media [20[. Related techniques which allow coherent focusing in scattering 
media are known from ultrasound j^lj and microwaves {22]. The frequency of those types 
of waves is low enough that electronic transducers can be used to time reverse waves, which 
redirects the waves towards their source. This technique has successfully helped to improve 
imaging resolution {23] and communication bandwidth 24]. 

In this Letter we report the first experimental demonstration of combined spatial and 
temporal control of light trajectories through opaque media via spatial wavefront shaping. 
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We apply our method to create an ultrashort pulse from the light transmitted through a 
strongly scattering medium. We can control the amount of time the optimized pulse stays 
in the sample and thereby select the path length of the light through the medium. The 
efficiency of our method is independent of the time delay, allowing for ultrashort pulses to 
be formed even from very long paths in the sample. 

The experimental realization can be summarized as follows (Figured]). An short-pulse 
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FIG. 1: Experimental setup. A beam from a mode-locked Titanium:Sapphire (Ti:Sa) laser is 
coupled into a heterodyne Mach-Zehnder interferometer. In the signal branch the beam reflects 
off a spatial light modulator (SLM) and is subsequently focused onto the sample. The light in the 
reference branch is frequency shifted by 40 kHz by two acousto-optical modulators (AOM) and 
then passes through a motorized delay line. The signal from the detector is filtered by a lock-in 
amplifier (LIA) and recorded by the Computer. 

light source illuminates a phase-only spatial light modulator (SLM), which can alter the 
phase of the light reflected from its surface. The SLM pixels are grouped into N indepen- 
dent segments each of which induces a controllable phase shift A$j, which can be considered 
constant over the bandwidth of the laser. The scattering sample is placed in the Fourier 
plane of the SLM. Both SLM and sample are embedded in the signal arm of a heterodyne 



Mach-Zehnder- type interferometer [251 ] . The heterodyne signal exactly corresponds to the 
cross-correlation of the of the forward scattered pulse with the reference pulse, which is 
delayed by a variable time delay r. The reference pulse is close to Fourier-limited, so that 
amplitude and phase of the transmitted pulse will not change significantly over the duration 
of the reference pulse. The heterodyne signal is then effectively an instantaneous measure- 
ment of the transmitted electric field. This signal serves as a feedback for an optimization 
algorithm, which programs the SLM. The time delay of the reference pulse is adjusted by an 
automated delay stage allowing optimization at a desired point in time r opt . An extended 
description of the setup including technical details has been included in the supplementary 



material. 
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FIG. 2: Optimized and random speckle pulses, (a) Amplitude of heterodyne signal of a non- 
optimized pulse as a function of time delay, averaged over 50 random speckle pulses, (b) Typical 
single random speckle pulse, (c)-(g) Amplitudes of single pulses after optimization at different 
time delays which are indicated by the dashed arrows. The optimization has been performed by 
dividing the SLM into 300 segments. The optimization generates strong, short pulses from diffuse 
light. The zero delay position is at the maximum amplitude with no sample. The plotted curves 
have been normalized to the maximum of the average non-optimized heterodyne signal (factor 1.53 
mV -1 ). 



The principle of the experiment can be described as follows. Light reflected from a single 
segment on the SLM is transmitted through the sample, giving rise to the field Ei(t) at the 
detector. Its phase can be modified by a time-independent phase shift A$j via the SLM. 
The total field scattered into the detector E out {t) is therefore given by the sum over all 
segments 

N 

i=i 

Multiple scattering allows us to assume that the contributions Ei(t) from the different 
segments at every single point in time t are uncorrelated random variables with Rayleigh 

rn. 

distributed amplitudes \Ei(t)\ and uniformly distributed phases [26]. For the non- 

optimized case, the resulting field E ont (t) can be viewed as the result of a random walk in 
the complex field plane. After the optimization, all contributions are in phase, adding up 
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constructively. The average amplitude enhancement is given by 17| 

(a) - 



I -^opt I ) rms 



( | -^rnd | ) rms 

= {l{N-l) + lfl^{\Nf/\ (2) 

The average intensity enhancement rj can be obtained by rj = a 2 . 

The non-optimized data was obtained by setting random phase values to the SLM seg- 
ments. The optimization algorithm adjusts the phase shifts A$j such that the amplitude of 
the heterodyne signal is maximized. We performed the optimization at 20 equidistant time 
delays r opt between -1.05 ps to +13.6 ps. For each r opt , the optimization was performed four 
times, with N = 12, 48, 192 and 300 segments, respectively, each time starting from a new 
random phase pattern. 

Our main result is displayed in Figure EJ showing the amplitudes of both the non- 
optimized (black lines) and the optimized (red lines) pulses for different time delays r opt 
and N = 300 segments on the SLM. The optimized amplitudes show sharp, distinct peaks 
with dramatically increased amplitudes at the desired time delay. We can control the amount 
of time the optimized pulses stay in the sample by the time delay r opt , and by that we con- 
trol the path length of the pulses through the sample. Note that the heterodyne signal is 
proportional to the field amplitude, the intensities exhibit even more pronounced optimized 
peaks. 
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FIG. 3: Enhancement a versus selected time delay r op t for different number of segments N on the 
spatial light modulator. 

The enhancement a versus time delay r opt is shown in Figure El Its magnitude, depending 
on the number of segments on the SLM, is constant from zero to several picoseconds time 
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delay. This result shows that our method works for short light paths as well as for light 
paths more than ten times longer than the sample thickness. 

For long time delays a continuous decrease of a is observed, which is related to the noise 
level of the experiment. We include a quantitative analysis of this effect in the supplemen- 
tary material. 

• Experiment 




FIG. 4: Average amplitude enhancement a (red dots) as a function of the square root of the 
number of segments N. The solid line is given by the expected a = 0.90 (^iV) 1 / 2 , without any 
free parameter. 

Figure H] shows the average enhancement in the constant regime in Figure E] versus N to- 
gether with the enhancement expected from theory (EqJ2J), (a) = cr(^N) 1 / 2 . The prefactor 
a = 0.90 corrects for the non-uniform illumination of the SLM surface with a truncated 
Gaussian beam, which effectively leads to a reduction of the number of used segments (see 
supplementary material). Our model matches the data well with no adjustable parameters. 

We have investigated the duration of optimized pulses in detail. Figure [5] shows the het- 
erodyne signal of three typical optimized pulses in a detailed scan of the time delay around 
the respectively chosen optimization time r opt . The average width of the optimized peaks 
(full width at half maximum) Ar opt = 190 ± 7fs, shows no dependence on time delay range 
between zero and 10 ps. The heterodyne signal is given by the amplitude of the optimized 
pulses, convoluted with the bandwidth-limited reference pulses. By deconvolution we cal- 
culated the field amplitude of the optimized pulses and determined their (intensity) pulse 
duration At opt . The optimized pulses have a duration At opt = 115 fs. For comparison, the 
input pulses have a transform-limited duration of At; n = 64 fs. This raises the question 
whether the lengthening of the transmitted pulse is caused by spectral narrowing or by re- 
maining fluctuations of the spectral phase. Numerical simulations, which are appended in 
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FIG. 5: Detailed cross-correlation scans around the optimized pulses at the time delays r opt = 1.1 ps 
(a), 4.8 ps (b) and 8.5 ps (c), together with a Gaussian fit (red lines). The amplitudes have 
been normalized to the maximum of the respective peak. The width (FWHM) of the peaks are 
Ar opt = 186 fs (a), 188 fs (b) and 206 fs (c). These values are typical in the range of time delays 
between zero and 10 ps, with an average width Ar op t = (190 ± 7)fs. 

the supplementary material, show that the answer is dependent on the number of segments 
N. For a low number of segments, the spectrum of the optimized pulse is random with a 
very short frequency correlation, but overall retains the shape and the width of the Gaussian 
input spectrum. The optimized spectral phase is not flat so that the pulses do not reach the 
bandwidth limit. For an increasing number of segments, the duration of the optimized pulse 
converges to its Fourier- limit. The spectrum has a smooth Gaussian shape, but a bandwidth 
narrower than the input spectrum. The method is capable of creating bandwidth-limited 
pulses, but since it is based on linear interferometry, the adaption to other pulse shapes is 
equally possible. 

The time-integrated intensity (energy) of the pulse with the highest peak depicted in Fig- 
ure |2] is 13.5 times higher than the energy of the average non-optimized pulse. In addition 
to the temporal optimization, overall more light is transmitted into the detected channel, 
demonstrating that the scattered light is controlled spatially and temporally. Our method 
exploits the mixing of spatial and temporal degrees of freedom by the random medium 



13J, to control the transmitted light in two spatial and one temporal dimension by only 



controlling spatial degrees of freedom on the two-dimensional SLM. On the one hand, the 
conversion of spatial degrees of freedom into temporal ones comes to the price of a speckle 
background, which on the other hand is easily outweighed by the enormous number of de- 
grees of freedom provided by state-of-the-art SLMs. The large number of controllable spatial 
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degrees of freedom is a great advantage over frequency domain pulse shaping techniques. 
Spatiotemporal control of the light field allows a far more generalized application of present 
coherent control schemes and marks a further step towards optical time reversal. 

In the experimental realization presented here, we optimized the pulse front using linear 
interferometry as feedback signal. The optimization of a non-linear response, like second- 
harmonic generation, will also lead to a comparably optimized pulse {27]. Using a nanopar- 



ticle with a non-linear emission response 



14| then enables the focusing of ultrashort pulses 



inside complex scattering media. We envision that our method can improve approaches for 
selective cell destruction in tissue [28]. In view of its potential for sharp focussing, it has 
potential for nanofabrication, nanosurgery and other micromanipulation techniques. 

Up to now we have not discussed the spatial extent of the optimized pulse. We use an 
aperture to select a single speckle spot in the Fourier plane of the sample for optimization, 
which corresponds to light transmitted into the forward direction. We know that transmit- 
ted fields in adjacent speckle spots are uncorrelated 29], from which we can conclude that 
the optimization is indeed limited to the selected area. An important future direction would 
be to investigate the spatial extent of the optimized pulse as a function of delay time. A 
combination with spatial scanning allows the measurement of the transmission matrix of the 
medium [30| in one temporal and two spatial dimensions. For Anderson-localizing samples 



311 ] . the size of the optimized speckle should be strongly time-dependent [32]. 

In conclusion, we have experimentally demonstrated that spatial wave front shaping of 
a pulse front incident on a strongly scattering sample gives spatial and temporal control 
over the scattered light. Our approach provides a new tool for fundamental studies of light 
propagation and has potential for applications in sensing, nano- and biophotonics. 
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